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The light curves of type I X-ray bursts (XRBs) result from energy released from the 
atmosphere of a neutron star when accreted hydrogen and helium ignite and burn ex- 
plosively via the rp-process. Since charged particle reaction rates are both density and 
very temperature dependent, a simulation model must provide accurate values of these 
variables to predict the reaction flow. This paper uses a self-consistent one-dimensional 
^ . model calculation with a constant accretion rate of M = 5 x 10 16 g/s (0.045M^.) and 
reports on the detailed rp-process reaction flow of a given burst. 

1. Introduction 

The rp-process in XRBs involves thousands of different reaction rates, which include 
(j), 7)-, (a, 7)-, and (a,p)-reactions, their respective inverse reactions, and electron cap- 
tures and /5 + -decays, and produces proton-rich and highly radioactive nuclei near or at 
the proton dripline pQ . Most of these reactions have not been measured, but with recent 
upgrades of existing experimental facilities and the construction of new facilities such as 
the proposed RIA and GSI-FAIR more reactions are becoming accesible to experimental- 
ists j2]. Yet radioactive ion beam experiments take a long time to prepare and execute, 
so in order to prioritize experiments it is important to accurately establish the reaction 
flow to determine which reactions are the most critical in determining the x-ray burst. 

With a neutron star surface gravity of g ~ 10 14 cm/ s 2 maintaining a hydrostatic pressure 
balance requires extreme densities and temperatures which are otherwise not obtainable 
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in terrestrial laboratories. Thus accreting neutron stars are also excellent nuclear labora- 
tories, because the emitted X-ray light curve depends indirectly on the nuclear reactions 
in different layers of the neutron star atmosphere j^j. 

The observed X-ray burst light curve obtains from the outbound heat transport of 
the nuclear burning front as it spreads around the star [J] . whence direct comparisons 
require multi-dimensional modeling which is still computationally unaffordable. Therefore 
comparisons between models and observations are restricted to gross burst features such 
as rise and decay times for bolo-metrically single peaked XRB light curves. Luckily the 
local nuclear reaction flow, which will be an important ingredient in next-generation 
computational models and a guide to future experiments, depends only on the local values 
of temperature, density and composition. Currently these values are best determined by 
self-consistent one-dimensional models. 

In this paper we use such a model which is detailed in jJJJ- It includes a 298 isotope 
network, diffusive heat transport and convection, and solves the full general relativistic 
equations on a co-moving metric with a conservative discretization. The computational 
domain is discretized into 103 zones with a pressure ranging from P = 5.1 x 10 20 erg- cm -3 
to P = 6.7 x 10 23 erg • cm -3 and it is bounded by a relativistically corrected radiative zero 
atmosphere and a realistic core boundary interface. 

2. The nuclear reaction flow 

This paper considers a neutron star accreting at a constant rate of M = 5 x 10 16 gs~ 1 .At 
this accretion rate ignition occurs in a layer where hydrogen has burned completely into 
helium j^J. Consequently the light curve obtains from the catastrophic detonation of 
a helium-rich layer and the subsequent conflagration of a layer of mixed H/He located 
above. Depending on the violence of the helium detonation and the burning time-scale 
and the strength of the mixed H/He burning a double peaked burst can obtain [7]. 

Ignition occurs in the ashes of the previous burst at the hydrogen fuel surface (see 
[HI), where the concentration of 4 He is maximal with mass fraction, X a = 0.57. The 
rapidly released heat spreads inwards and outwards demarcating two burning regions: 
One comprising helium burning where no hydrogen is present and the other comprising 
hydrogen burning in the form of the rp-process. 

2.1. Helium flash 

Unlike the rp-process, helium burning does not depend on temperature-independent 
/3-decays, whence the helium burning region rapidly adjusts within the tenths of a second 
it takes to reach peak temperature of T = 1.4 x 10 9 K. As the runaway exhausts 4 He, 
the triple-alpha slows down markedly, because the helium burning triple-alpha process 
is very sensitive to the temperature and the 4 He density. Other reactions include (ct,p)- 
reactions on previously generated stable ashes between 23 Na and 34 S, where the released 
protons capture on the alpha-chain nuclei 12 S, 16 S, 20 Ne, 24 Mg, 28 Si, 32 S, and 36 Ar until 
the Coulomb barrier prohibits further (a,p)-reactions on the resulting compound nuclei. 

The rapid temperature rise obtained from the helium detonation establishes a slightly 
super-adiabatic temperature gradient which causes a rapid rise in surface luminosity and 
maintains a convective region for 0.9 seconds This region extends all the way to the upper 
boundary of our computational domain and mixes its composition thoroughly. Since the 
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helium burns in about a tenth of a second, the subsequent energy release of the XRB 
obtains from the rp-process in the mixed layer of H/He above. 

2.2. Mixed H/He burning 

Prior to the burst the composition of this region ranges from freshly accreted material 
to processed material of previous bursts which has been advected down from the surface 
while burning steadily in the hot CNO-cycle. The heat of the convective bubble of the 
helium detonation rapidly mixes the entire region which consequently assumes the same 
composition. Therefore the end-product of the burning depends on the temperature 
attained in a given layer. 

The upcoming heat wave triggers the triple-alpha reaction and the 14 0(a,p) 17 F and 
15 0(a,7) 19 Ne breakout reactions. The former establishes the hot-CNO bi-cycle of PQ, 
which breakout depends on the 18 Ne(a,p) 21 Na, whereas the latter establishes a flow to 
21 Mg via 19 Ne(p, 7) 20 Na(jo, 7) 21 Mg, where it is blocked by photo-disintegration, because 
of the 21 Mg(p, 7)(7,p) 22 Al-equilibrium. Therefore the flow proceeds via 21 Mg(/5 + , Tu% = 
0.122s) 21 Na(p, 7) 22 Mg, where it branches going into either 22 Mg(/3 + , Ty 2 = 3.86s) 
22 Na(p, 7 ) 23 Mg(p, 7 ) 24 A1 or 22 Mg(p, 7 ) 23 Al(p, 7 ) 24 Si(/5+, T 1/2 = 0.102s) 24 Al. 24 A1 
captures a proton to 25 Si, which is again in blocked by photo-disintegration from the 
25 Si(jo, 7)(7,p) 26 P reaction, where a fast additional proton capture may create the short- 
lived 27 S, which decays to 27 P. Following the 25 Si(/3 + , Ti/ 2 = 0.634s) decay, the reaction 
path proceeds via 25 Al(p, 7) 26 Si (T1/2 = 2.23s), which is too long-lived for its decay to 
be relevant, so an additional proton capture leads to 27 P. Here the flow branches into 
either 27 P(/?+,T 1/2 = 0.260s) 27 Si(p, 7 ) 28 P or 27 P(p,7) 28 S(/3+,T 1/2 = 0.125s) 28 P. This is 
followed by 28 P(p,7) 29 S(/3+,T 1/2 = 0.187s) 29 P(p,7) 30 S(/3+,T 1/2 = 1.18s) where both of 
the (p, 7)-reactions on the sulfur nuclei have low Q- values, whence the chlorine compound 
nuclei are immediately photo-disintegrated. 

The time-scale of this flow depends on the arithmetic sum of the /3 + -decay half lives 
of the waiting points in the reaction path 9\. These waiting points are circumventable 
by the the (a,p)-process which makes the flow much faster pQ. (a,p)-reactions, which 
are very temperature-dependent, occur on the following isotopes: 18 Ne, 21 Mg, 22 Mg, 24 Si, 
25 Si, and 26 Si. At this point it becomes difficult for the alpha-particle to penetrate the 
Coulomb barrier of the target, but depending on the exact reaction rate of (a,p)-reactions 
on 28 S, 29 S, 30 S, 32 Ar, 33 Ar, and 34 Ar it may be possible to extend the flow further. Details 
of these reactions and their effect on the observable light curve are discussed in 

Following the flow from 30 S, it branches into either 30 S(o;,p) 33 C1 or 30 S(p, 7) (7,2?) 
31 Cl(/3+,T 1/2 = 0.200s) 31 S(jo,7) 32 Cl(/?+, T 1/2 = 0.298s) 32 S(p, 7 ) 33 CI or 30 S (T 1/2 = 1.18s) 
30 P(p, 7 ) 31 S(p,7) 32 Cl(/3+,T 1/2 = 0.298s) 32 S(p,7) 33 C1. From this point the flow con- 
tinues with 33 Cl(p, 7 ) 34 Ar(/3+,T 1/2 = 0.844s) 34 Cl(p, 7 ) 35 Ar(p, 7 ) 36 K(/3+,T 1/2 = 0.342s) 
36 Ar(jo, 7) 37 K(jo, 7) 38 Ca. Following its /3-decay and subsequent proton-capture the flow 
breaks into the p/-shell nuclei with 39 Ca(p, 7X7, p) 40 Sc(p,7) 41 Ti(/3+, T 1/2 = 0.080s) 
41 Sc(jo, 7) 42 Ti. Here the flow will branch again with either 42 Ti(/? + , 7\/ 2 = 0.199s) 
42 Sc(p, 7 ) 43 Ti(p, 7 ) 44 V or 42 Ti(p, 7 )( 7 ,p) 43 V(p, a) 44 Cr(/3+, T 1/2 = 0.053s) 44 V, where 
they connect and continue with 44 V(p, 7) 45 Cr(/3 + , T x / 2 = 0.050s) 45 V(p, 7) 46 Cr, which 
decays to 46 Ti via 46 V. 

At the lower temperatures near the surface the rp-process effectively stops here, but at 
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high temperature, the flow moves on via proton captures on these three A = 46 isotopes, 
which comprise the first bottlenecks to the heavier iron-group nuclei. The next bottleneck 
is the branching point 49 Mn which may either /3-decay or proton capture. In the latter 
case the reaction flow out of 50 Fe is limited by photo-disintegration, so the main flow 
proceeds with two proton-captures on 51 Fe and ends up on 53 Ni. The Z = 28 Ni-isotopes 
are well-bound, whence the proton-capture Q— values are low. Depending on the Q- 
values, the flow can proceed via 53 Ni(/3+, T 1/2 = 0.045s) 53 Co(p,7) 54 Ni(/3+, T 1/2 = 0.140s) 
54 Co(p, 7) 55 Ni(p, j)(j,p) 56 Cu(p, 7) 57 Zn otherwise the flow is stopped due to the long 
half-life of 56 Ni. An equivalent issue arises when the flow reaches 58 Zn, 59 Zn, and 60 Zn 
following 56 Ni(p, 7) 57 Cu(p, 7) 58 Zn and consecutive proton captures and /3-decays. The 
low proton-capture Q-values on 63 Ge, 64 Ge, 67 Se, 68 Se, 72 Kr, and 76 Sr cause immediate 
photo-disintegration of the compound nuclei, so the flow and the resulting energy gener- 
ation and burst decay luminosity is regulated by their comparably slow /3-decays. The 
maximum temperature is limited by the self-consistently determined hydrostatic pressure 
of the burning region and compared to the half-lives of the waiting points, the heat is 
quickly transported to cooler adjacent regions. Therefore the flow does not extend fur- 
ther except for trace amounts leaving an average composition of A ~ 63 and Z ~ 31 thus 
corroborating the results of |10j . 

3. Conclusion 

Matching observational luminosity curves to computational models require detailed 
multi-D models, which are still computationally unavailable. However, gross luminosity 
features and the nuclear physics of the burst can be determined from ID models. The 
burst rise time is determined by the characteristic time-scales of the burning regions. 
For mixed H/He burning, this is influenced by experimentally unknown reactions such as 
30 S(a;,p) 33 Cl and 34 Ar(a,p) 37 K. The surface composition following the burst depends on 
the proton-capture bottlenecks at A = 46 and 49 Mn. Finally the energy generation and 
composition of the burst ashes depend on the decays of the A > 64 waiting points whose 
distribution depends on the nuclear masses in that region. 
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